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TEMPORAL  EVOLUTION  OF  WHISTLER  GROWTH 
IN  A  COLD  PLASMA  INJECTION  EXPERIMENT 

Introduction 

Energetic  electron  precipitation  has  been  of  considerable  interest  for 
over  two  decades.  Bursts  of  precipitation  were  observed  by  balloon  x-ray 
bremsstrahlung  measurements  as  early  as  1963  [Winckler  et  al.,  1963; 
Anderson  and  Milton,  1964]  and  in-situ  measurements  of  greater  than  40-keV 
electron  fluxes  were  made  by  Injun  3  [O'Brien,  1964]  .  The  importance  of 
loss  cone  electromagnetic  cyclotron  (EMC)  instabilities  in  the 
magnetosphere  and  the  relevance  of  such  instabilities  to  particle 
precipitation  processes  was  first  pointed  out  by  Cornwall  (1965),  and 
independently,  by  Obayashi  (1965).  An  attempt  to  formulate  a  theory  of 
particle  trapping  and  precipitation  under  the  influence  of  EMC 
instabilities  was  made  by  Kennel  and  Petschek  (1966).  Kennel  and  Petschek 
studied  the  wave  growth  and  loss  rates  of  whistler,  ion-cyclotron  and 
magnetosonic  waves  in  a  finite  plasma  (i.e.,  magnetosphere)  and  derived  a 
limit  on  the  stably  trapped  particle  fluxes  in  the  radiation  belts. 
Subsequently,  more  rigorous  and  self-consistent  studies  of  these  processes 
have  been  conducted,  with  the  result  that  some  of  the  conclusions  of  Kennel 
and  Petschek  must  be  modified  (Etcheto  et  al . ,  1973).  Cocke  and  Cornwall 
(1967)  suggested  that  cold  plasma  played  an  Important  role  in  controlling 
the  wave  particle  interactions  in  the  radiation  belts  within  the 
plasmasphere .  Later  Brice  [1969]  pointed  out  that  cold  plasma  injected 
from  the  ionosphere  into  the  magnetosphere  could  produce  enhanced  whistler- 
mode  turbulence  and  associated  precipitation  of  energetic  electrons,  while 
Cornwall  et  al.  [1970]  recognized  a  parallel  process  for  the  ion  cyclotron 
waves  and  ion  precipitation. 

Brice  [1970],  and  3rice  and  Lucas  [1971],  have  suggested  that  a 

substantial  increase  in  the  energetic  electron  precipitation  could  be 
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achieved  by  injection  of  very  modest  amounts  of  cold  plasma  into  the 
radiation  belts.  A  hypothetical  experiment  in  which  VLF  EMC  noise  and 
proton  precipitation  could  be  enhanced  by  lithium  injection  was  analyzed  by 
Cornwall  (1974).  Similarly,  enhancement  of  VLF  EMC  noise  and  electron 
precipitation  by  barium  injection  was  studied  by  Liemohn  (1974). 

Under  certain  conditions  the  injection  can  introduce  a  new  limit  to 
the  number  of  stably  trapped  particles  [Kennel  and  Petschek,  1966] . 
Electrons  with  energies  below  a  certain  threshold  energy  Ec  will  continue 
to  be  trapped  and  those  with  energies  higher  than  Ec  may  be  precipitated 
provided  conditions  regarding  their  degree  of  anisotropy  and  relative 
abundance  are  satisfied.  For  ambient  cold  plasma  densities  the  threshold 
energy  Ec^  is  greater  than  the  thermal  energy  of  the  hot  particles,  and  the 
energetic  electron  density  can  increase  to  a  large  value.  With  the 

injection  of  additional  cold  plasma,  the  threshold  energy  is  reduced  from 
E£^  to  Ec2  and  electrons  with  energies  between  Ecl  and  E^  create  growing 
waves,  are  pitch  angle  scattered  and  precipitate  into  the  atmosphere.  This 
precipitation  may  produce  an  artificial  aurora  which  is  observable  from  the 
ground . 

Cuperman  and  Landau  [1976]  have  studied  the  theory  of  the 

electromagnetic  electron  cyclotron  (whistler)  instability  produced  by 
addition  of  cold  plasma  to  an  infinite  uniform  anisotropic  plasma.  Their 
formalism  disregards  the  finite  size  of  the  region  of  cold  plasma 
enhancement  in  the  magnetosphere  and  uses  a  stationary  model  for  the  cold 


plasma  density.  We 

extend  th 

e  Cuperman  and  Landau 

model  to  study  a  t 

varying  cold  plasma 

density 

injection  nc(t): 

and 

take  into  account 

finite  sire  of  the 

flux  tube 

by  studying  the 

wave 

gain  as  a  function 

both  the  spatial  and  the  temporal  growth  rates  of  the  waves. 


The  purpose  of  this  brief  report  is  to  make  semi -quantitative 
estimates  of  the  magnitude  and  the  duration  of  the  enhancement  of  VLF  EMC 
noise  and  electron  precipitation  effects  which  might  be  produced  by  a  cold 
lithium  release  experiment  in  the  parameter  range  accessible  to  AMPTE. 

Model 

We  assume  a  point  injection  of  neutral  lithium  gas.  The  neutral 

lithium  has  a  Maxwellian  distribution  and  is  allowed  to  expand  radially 

into  the  vacuum.  Solar  radiation  ionizes  the  expanding  lithium  gas  with  a 

time  scale  of  3000  seconds.  After  ionization  the  cold  plasma  is  frozen  to 

the  geomagnetic  field  and  its  volume  increases  in  time  as  a  geomagnetic 

field-aligned  cylinder  (see  fig.  1).  The  cylinder  has  a  radius  Rq  and  its 

2 

length  is  a  function  of  time  L(t),  so  that  its  volume  is  irR^L(t). 

To  estimate  the  number  of  cold  lithium  ions  in  this  field  aligned 
cylinder  we  also  modeled  the  evolution  of  the  neutral  lithium  injection 
cloud.  We  assume  a  release  of  100  Kg  of  neutral  lithium  and  a  yield  of  6% 
so  that  the  total  number  of  vaporized  lithium  atoms  is  ~  1.2  x  10  .  The 
initial  distribution  of  neutral  atoms  injected  at  a  point  R^  has  a 
distribution 


f0(v,  R,  0J 


°  377  5^R_RJ  exp  v2/vo2) 


with  thermal  velocity  Vq  ~  4.5  km/sec. 


(1) 


Integrating  over  the  initial  distribution  gives  the  neutral  lithium  density 
as  a  function  of  radius  from  the  release  point  and  time 
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(2) 


where  x^  is  the  ionization  time  for  neutral  lithium  atoms.  Initially  we 
are  interested  only  in  the  total  number  of  ions  in  the  flux  tube  and  not 
their  position  in  the  tube,  so  we  have  neglected  transport  of  ionized 
lithium.  The  rate  at  which  lithium  ions  are  created  is  given  by 


d_ 

dt 


n1(R, 


t) 


i-  n0(R,  t). 


(3) 


Assuming  a  finite  extent  of  the  whistler  wave  packet  across  the  magnetic 
field  of  Rq  ~  0(A^)(~  50  kms)  we  integrate  (3)  from  0  to  Rq  and  in  time 

from  0  to  •  to  estimate  the  total  number  of  lithium  ions  deposited  within 

23  _3 

a  sphere  of  radius  Rq  (N^  ~  5  x  10  cm  )  which  becomes  a  field  aligned 

flux  tube  since  the  ions  are  trapped  by  the  magnetic  field.  The  actual 

rt/xi 

rate  of  appearance  of  the  lithium  ions  in  the  flux  tube  is  N^fl  -  e  ). 

Additional  ions  are  created  outside  the  flux  tube  which  have  much  lower 
density  and  are  not  accounted  for  in  this  estimate.  Estimating  the  effects 
of  ion  transport  we  find  that  the  cylindrical  length  L(t)  increases 
linearly  in  time  at  a  velocity  Vq  in  both  directions  along  the  flux  tube. 
Thus  the  time  dependent  model  for  the  cold  ion  density  is  given  by 


nc(t) 


- - - r  (1  -  e 

t2VH*R0  i 


-t/x. 


8  , 


12.8  x  10  (  Rq  /Lq  j  [ 


,  1  ~  e 


-T 


-j  per  meter’ 


(5) 
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where  Lq  *  2VqT^  and  T  *  t/t^.  This  density  is  somewhat  higher  than  the 

average  flux  tube  density  because  the  neutral  and  ionized  lithium  is  peaked 

at  the  center  of  the  flux  tube.  Figure  (2)  shows  the  decay  of  the  cold 

lithium  density  as  a  function  of  the  normalized  time  T.  However,  the 

ionized  lithium  also  is  magnetically  trapped  along  the  field  lines  by  the 

mirror  force.  The  condition  for  the  low-  energy  lithium  ions  to  be  mirrored 

above  the  ionosphere  after  being  released  at  the  equator  in  the  altitude 

1  2 

range  of  5  to  8  rg  is  y  mv^  £0*02  ev  (see  Appendix).  Since  the  thermal 
energy  is  much  greater  than  0.02  ev  the  majority  of  the  lithium  ions  will 
remain  trapped  at  high  altitudes  in  the  magnetosphere  and  will  be  unable  to 
leak  into  the  ionosphere.  This  will  ensure  a  saturation  of  the  cold 
lithium  ions  for  large  T  at  a  value  larger  than  the  lowest  value  shown  in 
figure  (2).  Based  on  the  presence  of  magnetic  trapping  we  expect  the  cold 
density  will  saturate  around  T  ~  2  at  a  value  of  around  1.0  cm~^. 

We  now  discuss  the  critical  resonance  energy  criterion  in  the  presence 
of  a  time  dependent  cold  plasma  density.  It  is  well  known  that  wave 
particle  interaction  can  occur  when  the  resonance  energy  ER  is  greater  than 


or  equal  to  the  threshold  energy  Ec .  ER  was  defined  by  Kennel  and  Petschek 


[1966]  as 


er  =  7  mVR  ’ 


where  VT 


(d  -  n 


is  the  resonance  velocity.  Taking  into  account  the 


dispersion  relation  for  the  wave,  Cornwall  [1972]  calculated  a  minimum 


energy  condition  for  resonance  given  by 


Using  the  dispersion  relation  as  given  in  Cuperman 


where  A£  =  (Qg/orl)  *. 
[1974]  we  can  write 


A 

c 


2 

a 


a(t)0  ’ 


(8) 


where  a  =  kp,  p  =  vt^/^e>  a  =  1  +  nc(t)/ny  is  the  energetic  electron 

2 

density  and  8  =  nwKT/ (B  / 8 tt)  ,  where  K  is  the  Boltzman  constant  and  the  sub¬ 
script  w  refers  to  the  warm  electron  component.  Thus  we  can  rewrite  (7)  as 


_ 1 _ 

a(t)8A(l  +  A)2 


(9) 


where  Et^  is  the  thermal  energy  of  the  energetic  electrons. 

In  the  absence  of  cold  plasma  a  =■  1  and  Ec  is  a  constant.  This  is 
indicated  in  Fig.  (3)  top  line.  Here  we  took  the  density  of  the  ambient 
electrons  to  be  0.5  per  cubic  centimeter,  their  thermal  energy  Eth  ~  2.0 

_3 

keV,  a  magnetic  field  of  10  gauss,  the  degree  of  anisotropy 
A  -  0,50  and  S  0.07  .  With  the  addition  of  cold  plasma  E£  is  sharply 
reduced  especially  at  early  times,  thereby  enabling  a  very  large  population 
of  energetic  electrons  to  interact  with  the  waves,  enhancing  the  whistler 
mode  turbulence,  and  thus  as  a  consequence  of  pitch  angle  scattering 
precipitate  electrons  into  the  atmosphere.  Although  Ec  does  increase 
somewhat  at  later  times  it  is  still  much  lower  than  in  the  absence  of  the 
injected  cold  plasma.  Further,  as  described  earlier,  magnetic  trapping  of 
the  cold  lithium  ions  will  not  allow  the  cold  plasma  density  to  fall  as 
rapidly  as  indicated  in  rig,  (1),  thereby  ensuring  a  low  Z,  for  a 
considerable  length  of  time. 
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Results 


In  the  previous  section  we  saw  that  the  introduction  of  the  cold 
plasma  considerably  lowers  the  threshold  energy  and  thus  enhances  the 
growth  of  whistler  noise  turbulence.  Since  the  flux  tube  where  the  cold 
plasma  is  enhancing  wave  particle  interaction  has  a  finite  length,  the 
temporal  growth  rate  alone  cannot  determine  the  extent  of  the  wave  growth 
because  the  waves  propagate  away  from  the  region  of  rapid  growth. 
Therefore  we  consider  the  spatial  growth  rates  and  finally  calculate  the 
net  gain  G,  of  the  wave  during  a  single  pass  through  the  flux  tube.  The 
condition  for  a  large  effect  is  GR  >1  where  R  is  the  reflection  coeffic'  t 
for  the  wave  at  the  ionosphere.  The  waves  with  In  G  »  1,  a  more  s  're 
condition,  will  experience  a  significant  gain  in  the  wave  amp  Je 
independent  of  R,  leading  to  strong  whistler  turbulence  and  considerable 
pitch  angle  scattering  of  the  electrons  interacting  with  the  waves.  The 
gain  G  of  the  wave  is  defined  as 


G 


(10) 


8 


y(ai,k)L(t) 

Vg(u,k) 


g 


L(t)  <  4re 
L(t)  >  4rg, 


(11) 


where  y(u),k)  is  the  temporal  growth  rate  of  the  whistler  mode  instability, 
V  =  3u)r/3k,  is  the  group  velocity  of  the  whistler  mode  and  ra  is  one  earth 
radius.  Since  the  whistler  mode  is  excited  primarily  in  the  equatorial 
region  where  the  magnetic  field  is  nearly  uniform  and  straight,  we  have 
limited  the  length  L(t)  in  (11)  to  4ra. 
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In  order  to  evaluate  y(w,k)  and  V^(w, k)  for  the  resonant  energy  we 
solve  the  expression  (32)  of  Cuperman  and  Landau  [1974] 


a  0)  — 

-  - — =  ,  ^  «  «o/ne, 

1  -  U) 


and  the  resonance  condition 


to  -  1  +  a  v  ,  v  =  -  , 

R  R  Vth 


simultaneously.  Eliminating  to  from  (12)  and  (13)  we  obtain  a  cubic 
equation  for  'a.'  given  by 


a3  +  a8a  +  —  =  0 . 


We  solve  (14)  using  a  numerical  root  finder  for  'a'.  Then  'a'  is  used  to 
evaluate  the  group  velocity  V  from  (12)  and  the  temporal  growth  rate 
y/&e  from  the  expression  (33)  of  Cuperman  and  Landau  [1974].  With  these 
quantities  g  is  evaluated  for  the  resonant  waves. 

The  observed  differential  flux  for  radiation  belt  electrons  is  not  a 
Maxwellian  for  energies  higher  than  30  keV  which  is  a  typical  Ecfor  ambient 
cold  plasma  distribution  between  5  and  8  r  .  The  observed  flux  is  given 
by  a  power  law  of  the  form 


,3  electrons/ cm-  STE8.  sec.  keV, 


which  corresponds  to  the  energy  distribution  function  given  by 


J 


(16) 


> 


I 


I 


» 


» 


f (E)  =  C  E*(1  +  a\ 
o 


We  have  considered  two  distributions  of  the  following  form  for  the 
energetic  electrons  corresponding  to  quiet  and  disturbed  magnetospher ic 
conditions  respectively; 


f 

w 


m 


n  ( 


wl2itkT 


,3/2 

J  exp 


+  a). 

> 


(17) 


where  o  is  either  3/2  or  3  and  C3/2  ~  4.72  x  10"*1,  C3  =  6.7  x  10^  (private 
communication,  D.  Williams,  1983).  We  note  that  the  power  law  applies  only 
for  particle  energies  in  excess  of  30  keV  which  is  a  typical  value  for 
E£.  Above  Ec  the  stable  flux  is  limited  as  shown  by  Kennel  and  Petschek 
[1966]  and  others.  Below  Ec  the  ambient  flux  is  not  limited  by  the  wave 
processes  and  can  therefore  be  larger.  For  ease  of  calculation  we  make  a 
conservative  estimate  of  the  flux  in  this  energy  range  by  extending  the 
power  law  to  lower  energies . 

Figure  (4)  is  a  plot  of  g(T)  which  is  In  G(T),  as  a  function  of  T 
for  a  =  3/2  for  characteristic  particle  energies  spanning  from  20  to  40  keV 
in  steps  of  5  keV.  We  chose  plasma  parameters  such  that  A  =  0.75, 

6  =  0.07,  E(.|1  =  2  keV  and  the  density  of  energetic  electrons  =  0.5  per 
cc.  The  gain  curves  rapidly  increase  until  the  injected  electrons  fill  a 
flux  tube  of  length  4rg  and  thereafter  show  a  very  mild  decrease  due  to  the 
decrease  in  the  temporal  growth  rate.  The  group  velocity  shows  very 
little  variation  with  time. 


i 


I 
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Figure  (5)  is  a  plot  similar  to  Figure  (4)  except  here  a  ■  3.  The 

overall  behavior  of  both  figures  (4)  and  (5)  are  similar. 

Figure  (6)  is  a  plot  of  g  against  time  for  a  typical  energetic 

particle  energy  of  30  keV  for  various  values  of  8*  It  shows  that  the  wave 
gain  is  quite  sensitive  to  the  value  of  0.  For  higher  8  the  gain 

increases.  Therefore  if  magnetic  fluctuations  are  taken  into  account  8 
will  fluctuate  and  with  it  the  wave  gain. 

Overall,  we  find  a  significant  wave  gain  for  a  single  pass  through  the 
flux  tube  especially  for  lower  energy  particles.  Nevertheless,  at  25  keV 
the  gain  factor  is  already  large  enough  to  ensure  strong  turbulence  and 

strong  pitch  angle  diffusion.  For  lower  energies  the  gain  factor  becomes 
too  large  for  the  linear  theory  to  be  meaningful;  this  indicates  that  non 
linear  processes  (i.e.  strong  pitch  angle  scattering)  are  important  at 
these  energies.  Figure  7  shows  a  gain  versus  time  plot  for  lower  energies 
(4  to  12  keV).  The  gain  g(T)  increases  to  a  very  large  number  as  the 
energy  decreases  until  the  particle  energy  is  around  8  to  10  keV,  while  for 
still  lower  energies  the  gain  is  reduced.  This  decrease  occurs  for 
particle  energy  near  Ec> 

Conclusion 

Using  a  simple  time  dependent  cold  plasma  density  model  we  have  shown 
that  a  cold  lithium  injection  in  the  AMPTE  parameter  range  can  give  rise  to 
whistler  mode  turbulence  with  significant  gain  to  the  wave  amplitude  in  a 
single  pass  through  the  flux  tube.  whistler  mode  growth  gives  rise  to 
pitch  angle  diffusion  of  the  energetic  electrons  which  then  precipitate. 
As  energetic  electrons  are  removed  by  precipitation  new  electrons  drift 
into  the  flux  tube  via  gradient  and  curvature  drift  motion.  An  equilibrium 
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can  be  established  in  which  newly  supplied  radiation  belt  electrons  in 
excess  of  the  new  stable  trapping  limit  are  continuously  precipitated.  The 
rate  of  power  input  into  the  ionosphere  can  thus  be  estimated  roughly  as 
follows: 

A 

Power  Input  to  Flux  Tube  eff  ,  _eff,T  side  /1QN 

- - 5 - r — — = — r  =  n  kT  V,  - -  .  (18) 

Area  of  Ionospheric  Foot  oj  drift  A^oot 

2 

In  (18),  the  ionospheric  crossection  of  the  flux  tube,  A  £oot~(irRQ)BM/Bj , 

(Bj  and  BM  are  the  magnetic  fields  in  the  ionosphere  and  the  magnetosphere) 

and  the  meridianal  crossection  of  the  flux  tube,  A  .  ,  ~R  L(t)  give  a 

side  o 

4 

typical  value  for  the  ratio  x  ^  •  Using  a  reasonable 

s  f  f  2  ^ 

distribution  for  the  energetic  electrons  we  estimate  n^  ~10  cm  .  Very 
large  wave  gain  occurs  for  particles  of  energy  10  keV  which  have 

pff  _o  5  _2  -1 

KT  ~  10  ergs  and  a  ^ t~5xl0  cm  sec  .  Using  these  numbers  in  (18) 

—  9  —1 

we  find  a  precipitated  power  of  2  ergs  cm  “  sec  .  This  value  of  the 
precipitated  power  is  large  enough  to  produce  a  visible  aurora.  The 
precipitation  would  be  expected  to  continue  producing  a  visible  aurora 
until  the  injected  cold  plasma  in  the  flux  tube  is  lost  or  destroyed  by 
magnetospheric  convection  or  other  processes. 
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nc  =  0.0,  0  =  0.07,  A  =  0.5 


Fig.  (3)  A  plot  of  the  ratio  of  the  critical  energy  to  the  background  warm 
electron  energy,  Ec,  against  T.  The  top  line  corresponds  to  nc  - 
0,  0  =  0.07  and  A  =  0.5.  The  lower  curves  are  for  A  =  0.5  and 
0.75,  with  0  =  0.0  7  and  ng(T) . 


•  • 
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Appendix 


Here  we  calculate  the  minimum  perpendicular  energy  for  lithium  ions  In 
the  magnetosphere  to  be  trapped  at  high  altitude  by  the  magnetic  mirror 
force.  The  total  energy  of  the  lithium  ions  is  given  by 


_  1  2  „  G  Mm 

E  “  2  °  VH  +  pB  ‘  "r 


1  0 

where  u  *  mv^/B,  and  G  is  the  gravitational  constant.  Requiring 
conservation  of  energy,  the  net  energy  in  the  magnetosphere  must  be  equal 
to  the  net  energy  in  the  ionosphere, 


"  EI’ 


In  order  to  be  magnetically  trapped  in  the  magnetosphere  the  ions  must 


have  v^j  =  0.  Thus 


1  2  s 
2  m  V1M  l 


5  Mj--?-) 

I  M 

"5 


Taking  6  ev  for  the  net  gravitational  energy,  0.5  ev  for  the  parallel 
thermal  energy  in  the  magnetosphere  and  1  -  Bj/B^  ~  287,  we  find, 


I  m  V1M  ~  0,02  eV> 


for  trapping  to  occur. 
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